HYDRODYNAMIC LIMIT FOR A ZERO-RANGE PROCESS IN 
THE SIERPINSKI GASKET 
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Abstract. We consider a system of random walks on graph approximations 
of the Sierpinski gasket, coupled by a zero-range interaction. We prove that 
the hydrodynamic limit of this system is given by a nonlinear heat equation 
on the Sierpinski gasekt. 



1. Introduction 

The Sierpinski gasekt is a fractal in constructed in the following way. Start 
with an equilateral triangle of side 1. Divide it into 4 equilateral triangles of side 1/2, 
and retire the central triangle. Repeat this procedure on each of the three remaining 
triangles. After n — 1 steps, we are left with 3" small triangles of side 1/2". Since 
this sequence of triangles is decreasing, and each element in the sequence is compact, 
there is a non-viod limiting set K. We call this set K the Sierpinski gasket. 

At each step of this construction, consider the boundary of the resulting set 
as a graph r„ in M^, with 3(3" + l)/2 vertices and 3"+^ bonds. Consider now 
a system of particles evolving on this graph. The particles are attempting jumps 
to neighboring sites at rates that depend only on the number of particles sharing 
the same site. This is the so-called zero-range process. Such a system has been 
extensively studied in the usual lattice Z'^ or its periodic version (see [S] and the 
references therein). The purpose of this article is to study the collective behavior of 
this system as the graph get finer and finer, approximating in this way the Sierpinski 
gasket K. More precisely, we are interested in the hydrodynamic limit of the model, 
that is, the macroscopic evolution of the density of particles. It turns out that the 
hydrodynamic limit for this model is given by a nonlinear heat equation of the 
form dtu = A</)(u), where A is the Laplacian defined in K and (f> depends on the 
particular form of the interaction between particles. In a private communication, we 
have obtained Gaussian fluctuations for the hydrodynamic limit, in the stationary 
situation. 

At first sight, this result does not appear surprising, since the hydrodynamic 
equation for the zero-range process in the usual lattice Z'^ is also given by teh same 
nonlinear heat equation. However, the Sierpinski gasket if is a fractal, and therefore 
it does not have a differentiable structure. By this reason, even the definition of the 
Laplacian A in if is subtle. In [1], Barlow and Perkins have defined a Brownian 
motion in K as the scaling limit of the simple random walks defined in r„. As a by- 
product, they have defined the Laplacian A in if as the generator of this Brownian 
motion. A purely analytical definition can be found in the book by Kigami [B]. One 
remarkable fact is that the scaling limit is subdiffusive, in the sense that the scaling 
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factor, which is equal to 5" in this case, is larger than the square of the mesh of 
the graph, which is 2". This fact leads to anomalous diffusion properties of the 
corresponding heat equation, and to the failure of the usual Gaussian estimates on 
the decay to equilibrium of solutions of the heat equation. 

Although very detailed information about the fundamental solutions of dtu — 
Au in K has been obtained [7J, up to our knowledge existence and uniqueness 
of the Cauchy problem for the hydrodynamic equation dtU — A(f)(u) have not 
been considered in the literature. We provide in this article the existence and 
uniqueness results for the Cauchy problem in K that we need to make sense of the 
hydrodynamic limit for the zero-range process. Although collateral to our work on 
the zero-range process, these results could be of independent interest. 

From the point of view of interacting particle systems, the zero-range process is 
an example of a process satisfying the gradient condition. Roughly speaking, the 
current of particles is the gradient of another local function (the interaction rate, 
in this case). Therefore, Pick's law is satisfied at a microscopic level. The classical 
method to deal with the hydrodynamic limit of gradient systems was introduced 
in |4] and it is based on the so-called one-block and two-blocks estimates. Unfortu- 
natedly, the two-blocks estimate does not seem to hold for the Sierpinski gasket. In 
fact, this two-blocks estimate is based on the moving particle lemma, which states 
that a particle can be moved from one site to another with a diffusive entropy cost. 
The Laplacian A satisfies a Poincare inequality in K, and therefore, the spectral 
gap for the associated particle system is expected to be of the right order. However, 
in the Sierpinski gasket there are hot spots, that is, sites of the graph that have 
to be visited in order to connect two different regions of the graph. In particular, 
it is not true that the best strategy in order to transport a particle from one site 
to another is just to follow the shortest path, in contrast with the situation in the 
integer lattice. 

Probably the simplest method to prove hydrodynamic limits of gradient systems 
is the i/_i-norm method, due to Chang and Yau [2] (see ^ for a more readable 
exposition). The advantage of this method is that only requires the one-block 
estimate; its main drawback is that it works only for diffusive systems without a 
drift term. The idea of this method is very simple: let u, v be two solutions of the 
hydrodynamic equation. Then, the time derivative of {u—v, {—A)^^(u—v)) is equal 
to —2{u — V, (t>{u) — 4>{v)) and in particular, it is decreasing. Therefore, if uq = vq, 
then ut = Vt for any t > 0. The point is that this relation holds also at a microscopic 
level. However, at the microscopic level and for non-transitive graphs, there is a 
correction term involving AQ(x, x), where G is the Green function associated to A. 
Por the Laplacian in subsets of M'', Q{x,x) may not be well defined, but anyway 
AG{x,x) can be defined, since the singularity of Q{x,x) at the diagonal is always 
of the same magnitude, and cancels, at least in a weak sense, when taking the 
Laplacian as the limit of averaged differences around x. This is not the situation in 
the Sierpinski gasket K. We will see that, despite the fact that the Green function 
is continuous in K x K, the function G{x) — Q(x,x) satisfies AQ(x) = -t-oo, in a 
sense to be precised later. This nonregularity of the Green function poses an extra 
difficulty to the proof of the hydrodynamic limit. Usually, the one-block estimate 
allows us to replace a local function of the number of particles by averages over small 
boxes, when averaged with respect to a continuous test function. In our case, we 
need to average it with respect to discrete approximations of AQ(x), with we have 
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seen that is not continuous. Therefore, the one-block estimate needs to be proved 
without averaging with respect to test functions. This has been accompHshed only 
recently f5]. The proof however, confines us to dimension d < 2. Fortunately, the 
Sierpinski gasket has HausdorfF dimension dn = log 3/ log 2 < 2, and the so-called 
local one-block estimate holds in our case. 

This paper is organized as follows. In Section[2] we define the Laplacian operator 
in K and we review some aspects of functional analysis on K that will be needed 
in the sequel. In particular, we define with some detail the Green function Q{x,y) 
associated to the Dirichlet Laplacian in K. Although most of the material of this 
section has been taken from [6^, [Tj, we have decided to include most of the proofs 
for the convenience of readers interested in interacting particle systems, not familiar 
with analysis on fractals. 

In Section [3] we define what we understand by a weak solution of the hydro- 
dynamic equation dtu = A(/>(w) and we prove existence and uniqueness of such 
solutions by considering a finite-difference numerical scheme to approximate those 
solutions. 

In Section |4] we introduce the zero-range process, the iJ_i-norm method and 
we prove the hydrodynamic limit for the zero-range process, relying on the one- 
block estimate and suitable properties of the Green function Q{x,y). In Section 
[5] we prove the one-block estimate and the lemmas needed for the derivation of 
the hydrodynamic limit. In the Appendix we study the behavior of A„C/(a;,a;), 
where A„ corresponds to the discrete approximation of A defined in r„. The 
results presented in the Appendix are not needed for the hydrodynamic limit. Ee 
have included them here to stress that the example of the Sierpinski gasket is the 
worst possible case for the derivation of the hydrodynamic limit with the _ff_i-norm 
method. We believe that the ideas presented here can be adapted to treat general 
non-homogeneous graphs, like random graphs and trees, percolation clusters and 
disordered lattices. 

We have organized the paper in such a way that Sections |3] and |4] are independent. 
Therefore, the readers interested in the zero-range process can take for granted the 
existence and uniqueness results for the hydrodynamic equation and jump directly 
from Section [2] to Section 31 

2. The Sierpinski gasket 

Let Op = (0,0), ai = (l/2,\/3/2), 02 = (1,0) be the vertices of an equilateral 
triangle of unit side in M^. Define ipi : M.^ ^ by taking 'pi{x) — {x + ai)/2, 
i = 0,1,2. The Sierpinski gasket K is defined as the unique non-empty compact 
subset K of M? such that 

j=0,l,2 

A constructive definition of K is the following. Define Vn C recursively by 
taking Vq = {00,01,02}, Vn+i Uiipi{Vn)- Define V* = U„K- Since Vn C Vn+i, 
it is not hard to see that K — cls{V*), the closure of V* under the usual topology 

of 

Consider Vn as the set of vertices of a non-oriented graph r„ = (i?„,F„), and 
define inductively the set of bonds of r„ by taking Eq = {(ooai), (0102), {0.20,0)} 
and 

En+i = {{ifi{x)(fi{y)); (xy) G En,i^ 0, 1,2}. 
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Wc say that r„ is the n-th discrete approximation of K. For x^y G Vn, we say 
that X y if {xy) G E^. We simply write x ~ y when there is no risk of confusion. 
We say that x and y are neighbors in that case. 

2.1. The Laplacian operator in K. Let m : ^ R be an arbitrary function. 
For each n > 0, we define 

£n{u,u) = {h/ir E {<y)-<^)f- 

Proposition 2.1. For each n > and each u -.V* ^'R, 

£n+l{u,u) > £n{u,u). 

Moreover, given n > 0, m : y* ^ M, there exists a unique function Un V* ^ M. 
such that 

£n{u, U) = £n+p{Un, Un) for all p > 

and Un{x) = u{x) for every x G Vn- 

Proof. Let us consider the following setup. For a = (ao,ai,a2)) P = {Po^ Pi^ 
define the function 

Let a e be fixed. A simple computation shows that 

inf I{a, /?) = 5/3 \{ao - ai)^ + (ai - a2)^ + (0:2 - cto)"^] , 

and the infimum is attained in a single point /3, with Pi = {2a — ai)/b, where 
a = aQ + ai + 0L2- The proof follows easily from this observation and a chaining 
argument. □ 

From the previous result, £{u, u) = lim„ u) is always well defined, although 

maybe infinite. Observe that #y„ = 3(3"+l)/2. In particular, lim„ #y„/3" = 3/2. 
This motivates the following definition. For each n > define the positive measure 
fin in K by 

Unidx) = ^ 51 ^=o{dx), 

xev„ 

where dx{dx) is the Dirac mass at x. It is not hard to see that Hn converges in the 
vague topology to a measure n in K that coincides with the Hausdorff measure in 
K. A simple summation by parts shows that for any u:V*^'M., 

£n{u,u) = - j u{x)AnU{x)lJ,n{dx), 

where A„ is the discrete Laplacian in Vn- 

Anu{x) = 5^ J2 Hy)-u{x)]. 

Proposition 2.2. There exists a universal constant c > such that for any u : 
— > R with £{u, u) < +00 we have 

luiy) — u(x)\ ,1/9 

x,yev' \y ^1 
x^y 

where a = log(5/3)/21og 2. 
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Proof. For each x,y G Vn, define 

\u{x) - 



Rn{x,y)= sup 



£„{u,u) 

By definition, \u{x) — u{y)\ < ^jRn^x, y)£n{u, m)^/^. Notice that for any con- 
stants a > 0, 6 e M we have £„(au + b,au + b) ~ a^£„(u, u). Therefore, 

Rn{x,yy^ ^ inf £„(u, u). 

u{v) = l 

Assume that x V- Considering the function u{z) = l(z = y), we see that 
Rix,y)-^ < 4(5/3)". By definition, R„{x,y)-^ > (5/3)". Therefore, 

\u{x) ~ u{y)\ < (3/5)"/2£:„(u,zi)i/2 < 2-""£:„(7.,7.)i/2, 

where a = log(5/3)/2 log 2. Remember that |x — y| — 2^" for x ^„ y. Therefore, 
we have proved the inequahty for x y when some n > 0. Using the triangle 
inequality, we can extend this relation to x,y € V* arbitrary. □ 

In particular, the previous proposition tells us that any function u : — > M with 
£{u,u) < +O0 is uniformly continuous. Therefore, u can be continuously extended 
in a unique way to the set K. From now on, we consider u as defined in K, and 
we will assume that u : K ^ R with £{u, u) < +oo is continuous. Notice that the 
points a^, i — 0, 1,2 are different from the other points in V* . In fact, the points 
Qi have only two neighbors in r„ while other points in Vn have 4 neighbors. It is 
natural to define Vq as the boundary of K. Define 

H^iK) = {u: K ^R; £{u, u) < +oo, u{a,) = Vi}. 

For u e H'^{K), define ||u||i = £{u,uY^^. Notice that £{u,u) = if and only 



if u is constant in K. Therefore, || • ||i is a norm. By Proposition 2.2 Hi{K) is 
closed under this norm. It is easy to see that {Hi{K), \\ ■ ||i) is a Hilbert space, 
with inner product given by the polarization identity 

£{u, v) = {£{u + v,u + v) — £{u — v,u — v)) /4. 

Denote by Co{K) the set of continuous functions u : K ^ M. with u{ai) = for 
every i. We define C^ilJ.) as the completion of Co(-R') under the norm 

iHio = { I u{xrf,{dx)Y^'. 

Proposition 2.3. The space H^{K) is dense in C'^ifJ.)- 

Proof. Notice that H^{K) C Co{K) C C^iK). It is enough to see that H^{K) is 



dense in Co{K). Take u G Co{K), and define m„ as in Proposition 2.1 Notice that 
given three numbers ao, ai, Qf2, Pi = {2a — ai) /5 is between the maximum and the 
minimum of a^. Therefore, the function Un satisfies a maximum principle, and for 
every x V* , 

\u{x) - Un{x)\ < \u{x) - u{y)\ + sup \u{x') - u{y')\ 

x' ,y' 

x'^nV' 

for some y E Vn. Therefore, 

sup \u{x) — Un{x)\ < 2 sup \u(x) — u{y)\^ 

xeV' x,y£V* 

\x-y\<2-" 
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which goes to as ti ^ oo. □ 

Denote by {u,v) the inner product in L^{^). Define in an analogous way the 
spaces C^{fj,n) and the inner products u), {u,v)n. Remember the identity 

£n{u,u) = -{u,AnU)n- 

By analogy with the case of the real line, we define the Dirichlet Laplacian 
A : D{A) C C^{fJ,) as the unbounded operator given by 

i) ^(A) = {u e H°{Ky,3c > with £{u,v) < c||w||o Vw G H°{K)}. 

ii) For u £ -D(A), Au — h if and only if £{u, v) — -~{h, v) for every v e Hq{K). 

Notice that this definition of the Laplacian in K is not constructive. At this 
point, even to find a single example of a function u G -D(A) is hard to achieve. 
Moreover, for a generic u € Co{K), the approximations Un of Proposition 2.1 are 
not in D{A). For this reason we will extend the definition of the Laplacian in the 
following way. For u € define the dual norm 

|2 

L — 

veH°{K) 



\u\\_i= sup \^2{u,v) — £{v,v)} . 



By Proposition 2.1 for any u E Hi{K) we have HitHoo < c||u||i. Therefore, 
the Friedrich's inequality \\u\\q < c||m||i holds, and in particular ||u||_i < +oo for 
every u G C^{K). We denote by H_i{K) the closure of C^{K) under this norm. We 
extend the definition of A to the operator (still denoted by) A : Hi{K) H^i{K) 
such that 

{Au, v) = £{u, v) for all tt, u G H'^. 

Notice that A is now well defined as an element of H^i{K) for every u G Hi, 
and A is an isometry from H^^K) to H_i{K). Since C^{n) is dense in H_i, we 
conclude that D{A) is dense in C^{n)- For any function u £ Co{K), we have 

AUn = ^ u{x)6x, 
x£V„ 

and therefore we are plenty of examples of functions u for which Au can be evaluated 
(in H^i, of course). Since C Co{K), the set A4'^{K) of Radon measures in K\Vo 
is contained in H^i(K), and the previous expression makes sense. We say that Un 
is the harmonic continuation of u\v„- Notice that Un converges to u in Hi{K), and 
therefore Au„ converges to Au in H^i{K). 

2.2. The carre du champ. The set K does not admit a differentiable structure, 
since it is clear that no neighborhood of K can be put in bijection with an open set 
of M'' in a differentiable way. Therefore, the notion of a gradient Vu for functions 
u : _ftr — > M seems hopeless. Moreover, since the dimension of K is not an integer, it 
is not clear how many components should Vu have. What is remarkable, is that the 
so-called carre du champ |Vup can be defined in a very simple way. We say that 
a set T (- K IS a, triangle if T = if U A{xq,xi,X2) for some triangle A{xq,xi,X2) 
with vertices mutually adjacents in Vn for some n > (that is, Xi Xj for i ^ j). 
In an equivalent way, T is a triangle if T is of the form Lpi^^ o • • • o ip^_^ [K) for some 
sequence {ii, . . . , j„} in {0, 1, 2}. Fix a function u e Hi{K). For each triangle T, 
we define 

(^M[n,«] = lim (5/3)" ^ (u(2/) - u(x))^ 



^ x,yeT 

V 



HYDRODYNAMICS ON A FRACTAL 



7 



By the proof of Proposition |2.H this sum is increasing and therefore the limit 
always exists. Moreover, the limit is always finite, since it is bounded by £(u,u). 
The set of triangles generates the Borel topology in K. It is also not hard to 
check the continuity at vacuum of the set-valued function fi[u,u] ■ Therefore, we 
conclude that ^j,[u,u] can be extended to a positive, finite measure in K. For two 
given functions u, v in Hi{K), we define the measure ^[u^v] by polarization: 

^[u,v\ '^{^^[u+v .u+v\ ^[u—v^u—v])- 

It has been shown that the measures are singular with respect to the 

HausdorfF measure /i for any u e Hi{K). However, it has been shown that these 
measures are not mutually singular in the following sense: there exists a measure 
fl va. K such that for any pair of functions u, v in Hi{K) we have 

for some function T{u, v) in £^(/2). In fact, the measure p. can be chosen as equal to 
+ /^[/i2-'i2] for suitable harmonic functions hi, h2- Therefore, we can define 
Vu • Vw = r(u, v) to get the identity 



£{u,v) = J Vu-Vvdfj,. 



2.3. Harmonic functions and the integration by parts. Usually, a function 
h : K R is said to be harmonic if Ah — 0. Notice that the only function h in 
Hi(K) for which Ah = is = and we only have defined Ah for functions in 
Hi{K). Therefore, we need a definition of what we mean by an harmonic function. 
A function /i : F* ^ M is said to be harmonic if A„/i(x) = for every a; € \ Vq 



and every n > 1. In that case, 8{h, h) — £o{h, h) and by Proposition 2.1 h can be 
uniquely extended to a continuous function h : K —t 'K. Notice that h is entirely 
determined by its values at the boundary Vq. 

We extend the definition of A as follows. For a continuous function u not nec- 
essarily in H^{K), we say that Au = u if there exists an harmonic function h such 
that u — h E H^{K) and A{u — h) — v. In this case we say that u e Hi{K). 

For a function u : — > M, we define the discrete Dirichlet Laplacian by 
A^u{x) = Anu{x) if a; ^ Vb and A^u{x) = if u G Vq. Define the normal 
derivatives d^^u by 

5>=(5/3r Yl 

We have the following (discrete) integration by parts formula: 

j=0,l,2 

In order to obtain an analogous formula for the Dirichlet Laplacian, for u G 
H^{K) we define 

d'u= {Au,h'), 

where h^ is the harmonic function with h^(aj) — Sij. Since A is symmetric, 
{u,Av) — {v,Au) for any pair of functions u,v & H^{K). It is straightforward 
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to check the identity 

(u, Au) - (w, Au) = ^ {u{ai)d'^v - v{ai)d''u} 

1=0,1,2 

for any two functions u,v Cz Hi{K). 

2.4. The Green function in K. Friedrich's inequality tells us that the Dirichlet 
Laplacian has a positive spectral gap in C^{^). In particular, for every u e 
the equation 

{—Aw = u ("2 11 

has a unique solution. Since the inclusion Hi{K) C C^{ij) is compact, we conclude 
that the operator (— A)~^ is compact. We put the — sign to emphasize that A is 
non-positive. In particular, there are an orthonormal basis {vi\i of C^{n) and a 
non-decreasing sequences {Xi}i of positive numbers such that — Awj = A^u^ for any 
i. The function w can be written in terms of the orthonormal basis {vi}i: 

i>i 

Formally, we can obtain w{x) by an integral formula: 

w{x) = J g{x,y)u{y)n{dy), where g{x,y) ^ ^ A,^^Ui(x)wi(y) 

i>l 

is the Green function associated to the Dirichlet Laplacian A in K. A simple 
computation shows that the sum defining Q{x,y) is convergent in (E) /i) if 

E,>i < We will give a constructive definition of g(x,y) that will allow 

us to prove finer properties of Q{x,y). Take the discrete Laplacian A„ in V„ and 
define the Green function Qn{x^y) at x^y S Vn as the solution of 



AnQn{x,y) 



0, x£ Vb, 



where 5{x, y) — 1 ii x — y and 5{x, y) — otherwise. Here the operator A„ acts on 
the first variable x. Notice that Gn{x, y) is non-negative and Qn{x, y) < Gn{x, x) for 
any x,y G Vn- We extend the definition of Qn{x,y) to K by taking the harmonic 
continuation of Qn given by Proposition |2.2[ A key observation is that for y ^Vn, 
Gn+i{.x,y) = Qn{x,y) for any x G Ki+i- In fact, for x e K+i with x ^ y, 
An+iGnix, y) = 0. For x = y, a, simple computation shows that Gn scales correctly, 
and therefore An+iGn{y,y) = 3""'""'^. The following propositions shows that it is 
sufficient to compute Gi{x,y) to obtain Gn(x,y) for every n > 1, x,y € V*: 

Proposition 2.4. For any i = 0,1, 2, 

Gn+i{(Pi{x),Lpt{y)) = -gn{x,y) + ^ h^{y)gi{(pi{x),ipi{a.i)). 

J=0,l,2 



The proof is simple; we refer to Section [A] for the argument. As a consequence 
of this relation and the previous discussion, Qn{x,y) does not really depend on 
n. Therefore, we define for x,y G V*, Q{x,y) = Qn{x,y), where n is such that 
x,y & Vn- For fixed y, we see that £iG{-,y),G{-,y)) = Q{y,y)- In particular, 



HYDRODYNAMICS ON A FRACTAL 



9 



G{',y) is uniformly continuous and can be uniquely extended to K. In this way we 
can not define Q{x,x) for x ^ V* . We would like to prove that in fact Q{x,y) is 
uniformly continuous in V* xV* . This is an immediate consequence of the following 
proposition: 

Proposition 2.5. There exists a constant c > such that Q(x,x) < c for every 
X e V*. 

Proof. Notice that ^^h'^ix) = 1 for every x ^ K. In fact, J^iJ^i corresponds to 
the harmonic function with h{ai) = 1 for i = 0, 1,2, which is identically constant. 
Taking x = y in Proposition |2.4| we see that 



G{ipi{x),(pi{x)) = -g{x,x) + ^ h^y)g{ipi{x),ai). 

Since every y e Vn+i \ Vn is equal to (pi{x) for some i and some x G Ki \ Ki-i, 
we see that 

sup g{y,y)<i/5 sup g{x,x)+ sup sup g{x,ipi{a-i)). 

Setting f3n = s'^Pxev„\v„-iGi^^^)^ ^^e that Pn+i < 3/5/?„ + c for some 
constant c independent of n. A simple computation shows that /3„ is bounded in 
n. □ 

Remember that for any function u G Hi{K), converges to Aw in H-i{K). 
Therefore, AC/(a:,y) = 5y{x) and w{x) = J Q{x,y)u{y)fj,{dy) is the solution of 
equation ( |2.1[ ), at least for functions u € Co (j^r). By the continuity of Q{x,y), we 
conclude that J Q{x,y)u{y)fi{dy) solves (2.1) for u G C^{lj) as well. 

3. The nonlinear heat equation in K 

Let (f) : M+ M4. be a smooth function. We will assume that there exists eo > 
such that eo < < ^ ^o^' ^^^^ ^ ^ ''^+- -f'i^ some T > 0. We want to study 

the Cauchy problem 

dtu = A(j){u) 

u(t,ai) =ai,i = 0,1,2 (3.1) 
u(0, •) =-uo(-)- 

More precisely, we want to obtain criteria for existence and uniqueness of solu- 
tions for this equation. Now we define what do we understand by a weak solution 
of ( [3A| . We say that u : [0,T] x K is a, weak solution of ( [3l) if: 

i) For almost every t G [0,r], u{t, •) G Hi{K) and 

\\u{t,-)\\ldt < +00. 

ii) For any function G : [0,T] H^{K), pointwise differentiable in t and 
strongly differentiable in H^i{K) as a function of [0,r], 

{uT,GT)-{u^,Ga)- r {{uudtGt) + {c^{ut),AGt)]dt^ J2 c^^^'Gt. 

1=0,1,2 

We will start by considering a finite-difference scheme that approximates equa- 
tion dO. 
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3.1. A discrete nonlinear equation. Take a function Uq : Vn — + [0, oo) such 
that UQ(aj) = ai. Let us define u"-{t,x) : [0,oo) x F„ — > [0,oo) as the solution of 
the following system of ordinary differential equations: 

= A„0(u"(i,a;)) for x G V° 

u"'{t, Qi) — tti for i — 0,1,2 
u"(0,a;) ^ u^{x). 

By the maximum principle and Peano's theorem, u^{t,x) is well defined for any 
t > 0. We will prove existence of solutions for equation ( |3.1| in a proper sense 
by taking limits of these approximated solutions u^{t, x). To avoid an overcharged 
notation, we will take ai = 0. Our arguments work for Ui arbitrary as well: just 
take into account the boundary terms when performing integrations by parts. Let 
us define the discrete norms 



3'^ 

x~„y 

Let us denote the function •) by u". It is not hard to see that ||u"||o,n is 
decreasing. In fact, 

|||<||^,„ = 2«,A„0«))„ 

= -2£„«,0«))<-2eo||<||l„- 
Integrating this inequality between t = and t — T, we see that 

,n||2 I r,^ j IL,n||2 J-^^IL, l|2 



|K||^,„ + 2eo/ M\\indt<\\uo\\i,,. (3.2) 
Let us fix some reference time T > 0. For a function u : [0,T] x K ^ M., define 

IIHII?= TlHI?, 

and denote by ° ri-^) the Hilbert space obtained as the closure of C ( [0 , T] , Hi (K)) 
with respect to this norm , where C([0, T], Hi{K)) denotes the space of continuous 
paths in Hi(K). Fix a continuous function uq : K —> [0,oo) with u{ai) — for 



i = 0,1,2. Consider u", the harmonic continuation of u" into K. By (3.2), we have 

sup||K|||?<^. 

In particular, there is a subsequence n' such that u" converges to some function 
ut S Hi^t{K), weakly with respect to the topology of Hi t{K). 



Theorem 3.1. The function Ut € Hi_t{K) is a weak solution of (3.1) 



Proof. Taking a second subsequence if necessary, we can assume that 4>(u^) con- 
verges to some function Wt as well, where (t>{u^) is the harmonic continuation of 
4>{Uf). At this point, we need to justify the identity Wt — 4i{ut). For the usual 
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Laplacian, defined on a bounded, open set J7 C M'*, the argument is the follow- 
ing. If w e H1{U), that is, if (u, — Am) < +00, then there exists a unique func- 
tion Vu : [/ ^ M'' such that (u, AG) = (Vw, VG) for any smooth function G. 
Since 4>{u) is a smooth function of u, <j){u) also belongs to Hi{U), and moreover 
V0(u) = (j)'(u)\Iu. The uniqueness of the weak gradient Vu would allow us to 
conclude that wt — 4>{ut)- 

Recall the definition of the carre du champ Vu. A simple Taylor expansion shows 
that 



r 

as expected. Therefore, we can appeal to the uniqueness of the representation 
/^[u,ti] — r(u, u)]! to conclude that Wt — 4>{ut). 

Take a function Gf G Hi t{K). Assume that Gt is of class in time. By 
hypothesis, 

0«(a;))AGt(x)Ai(da;)ds / /" 0(ut(a;))AGt(a;)At(da;)di. 

Jk Jo Jk 

Performing an integration by parts, we see that the left-hand side of the previous 
expression is equal to 



T r rT 

dt 



Gt{x)^<l,{u^{x))^i{dx)ds^ / {^U^,Gt)ndt 

JK Jo 



{u^, Gt)„ - {uo, Go)„ - / «, 9tGt)„dt. 

Jo 



Remember that ||it^||o,n is also bounded by ||uol|oo- In particular, choosing a 
further subsequence if necessary, we can assume that converges weakly to ut 
in C^in). Therefore, {u'^,GT)n converges to {ut^Gx)- By Friedrich's inequality, 
weak convergence in Hi^t{K) is stronger than weak convergence in C^{fj,{dx) x dt). 
Therefore, we can pass to the limit in each of the terms on the right-hand side of 
the previous expression. We have therefore proved that 

{ut,Gt)~{uo,Go)- [ {{ut,dtGt) + {Hut),AGt}}dt = (3.3) 
Jo 

for any function Gf smooth enough, which proves the theorem. □ 



Theorem 3.2. The equation (3.1) has at most one weak solution. 

Proof. The heuristic argument is very simple. The idea is to formalize the following 
formal computation. Take two weak solutions ut, Vt of (3.1). Then, 

- vt\\\ = 2{{-A)-\ut - Vt),A{<l>M - <l>{vt))) ^3^^^ 

= -2{ut-vt,^iut)-(l){vt)) < 0. 

Therefore, if uq = Vg, we conclude that Ut = vt for any t > 0. Of course this 
heuristic computation needs to be justified. By (3.3 1, we have 

{ut-vt,Gt}^ I {{ut-vudtGt) + {^{ut)-<ty{vt)AGt)]dt. 
Jo 
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Let US take Gt — (—A) ^{ut — vt) in the previous expression. That is, define 
Gt{x) = / g{x,y){ut{y) -vt{x))fi{dy). 



IK 

Putting this into the previous formula, we obtain immediately 
||uT - -J^tII-i = -2 / {ut ~ vt,(l){ut) - (j){vt))dt, 



which is just the integral version of ( 3.4 1 . But we still need to justify that Gt can be 
taken as a test function. Since Ut and Vt are in Hi{K) for almost every t G [0,T], 
for fixed t the function Gt is regular enough. The problem is that Gt maybe is not 
differentiable in t. This problem solves easily by taking an approximation of the 
identity 75(t) with support in [0, ^] and defining 



Jo 



Gt+sds. 

Now Gf is differentiable, so it is an admissible test function. Taking 6 ^ we 
obtain the desired result. □ 



4. Hydrodynamic limit for the zero-range process 

4.1. The zero-range process. Let g : Nq = {0,1,- ••} [0,oo) be a function 
with g{0) = 0. The zero-range process in Vn with interaction rate g{-) is defined as 
the continuous-time Markov chain ^t in il„ = Nq" and generated by the operator 

Llr= E E5(e(^))[/(r'^)-/(0], 

where ^ is a generic element of il„, / : ri„ M and ^^'^ is given by 



Notice that the number of particles in this process is preserved by the dynamics. 
Therefore, for any fixed initial configuration, the state space is finite, and the 
previous process is well defined. This process has a family of invariant measures 
which we describe as follows. Define g{n)l — g{l) ■ ■ ■ g{n), g{0)l = 0. Assume that 




0* = I lim sup y^g{ny.\ 



is non-zero. This is the fact if, for example, inf„>„(, g(n) > for some uq. For any 
(j) < <j)* , define the uniform product measure Dfj, in ri„ by 

where Z{(j)) is the normalization constant. Notice that due to the fact that (/><(/>*, 
the normalization constant Z{(j)) is finite. It is not hard to see that the measure 
is invariant under the evolution of ^j. Observe that (j) — J g{^{x))i'^{d£,). 
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Define the number of particles per site by — J (^{x)D^{d^). Notice that the 
application i— > p(0) is strictly increasing, with p(0) = 0. Therefore, (p i— > p(0) is 
a bijection between [0,(f>*) and [0,p*), where 

p* = lim p(d>). 

Denote by p i— > (f>{p) the inverse mapping of p{(t>). Since the number of particles 
per site is a more natural quantity than (j), we define Vp — iy^(p) for p € [0,p*). 

Now we will introduce a Dirichlet-type boundary condition into this process. Fix 
some numbers G [0, p*), i = 0, 1, 2. Define the boundary operators by 

where the configurations £, i Sy are given by 



CW±1, z = y 



The zero-range process in Vn with boundary conditions {a^ji is then defined as 

v° 

the continuous-time Markov process in il„ = Nq" and generated by the operator 

w(o=E E5(e(:^^))[/r^")-/(0]+ E ^-/(^)- 

xeV° y£V° 4=0,1,2 

The difference between this process and the zero-range process defined previously 
is easy to understand. Inside Vn (that is, in V^), the dynamics is the same. Particles 
are coming from the boundary sites with intensity (/'(ai), which corresponds to 
have a density of particles ai at a^. Particles are also annihilated when they jump 
into the boundary sites a^, in order to keep the density of particles at fixed. 

Notice now that the number of particles is not longer fixed, since particles are 
coming in and out at the boundary sites. In order to have a process with amenable 
properties, we will impose some technical conditions on the interaction rate g{-) 
(see [5]). We say that g{-) satisfies (SG) condition if 

i) supfc \g{k + g{k)\ < +oo, 

ii) There exist fcg > and aq > such that g(k + l) — g{k) > a2 for any I > fcg. 

Notice that in this case p* — +oo. This condition guarantees the existence of 
exponential moments for the occupation variables ^{x) under the invariant measures 
Vp. We say that g{-) satisfies (C) condition if g{k + 1) > g{k) for any k. In this 
case, there exists a constant > such that / exp{9o^{x)}di'p < +oo for any 
p < p* . We will assume throughout this article that g{-) satisfies (SG) or (C). 

Condition (SG) guarantees the existence of a uniform spectral gap, which states 
that the magnitude of the first non-null eigenvalue of L^^ with respect to Vp is 
bounded below by a constant that does not depend on the density p. We expect 
this constant to be of order 5^". The non- validity of the moving particle lemma, 
discussed in the Introduction, prevent us to obtain such a bound. Notice, however, 
that a simple computation shows that there exists a constant cq, independent of n 
and p, such that the spectral with respect to Vp gap is bounded below by co6~". 

Condition (C) implies that the zero-range process is attractive, which means 
that, given two initial configurations ^, ^' with ^(x) < S,'{x) for any x, there exists 
a joint process {^t,it) such that is a zero-range process starting from ^, is a 
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zero-range process starting from ^' and ^t{x) < Cti^) for any x G Vn and any t > 0. 
This property allows us to obtain moment bounds for ^t(x) in terms of the invariant 
measures Vp. 

Now a simple path argument shows that there is exactly one invariant measure 
for the evolution of ^t. Remarkably, this invariant measure is still of product form. 
Consider the harmonic function h with h{ai) — ai. It is not hard to see that the 
non-uniform product measure vy^ defined by 

is invariant and ergodic for the evolution of ^j. 



4.2. The iJ_i-norm method: heuristics. Probably the simplest method to 
prove hydrodynamic limits for particle systems of gradient type is the i7_i-norm 
method introduced by Chang and Yau [5] (see |3] for a more comprehensible ap- 
proach) . The main drawback of this method is that only works for strictly diffusive 
systems. The other alternatives are the so-called entropy method |3] and relative 
entropy method [TO] . As we discussed in the Introduction, The entropy method 
requires a path lemma that roughly states that we can move a particle from one 
site to another paying a diffusive cost. This is not true for the Sierpinski gasket, 
due to the presence of hot spots: points that connects two huge parts of the graph 
that can not be avoided in order to move a particle from one of these parts to the 
other. For example, if we want to transport a particle from a site in Lp{){K) to other 
site in ifi{K), the particle has to pass by point c/?o(ai), or by points (^0(02), ¥'1(02). 

The second alternative requires smoothness of the solutions of the hydrodynamic 
equation. Of course, since we do not have a differentiable structure in K , we do 
not expect the solutions of the hydrodynamic equation to be smooth. 

The iJ_ 1 method is based on the heuristic argument leading to uniqueness of the 



hydrodynamic equation (3.11. Remember that the idea was to prove that, for two 
solutions ut, vt of (3.1), ||Mt — i"*!!?.! is decreasing in time. The main point is that 
this inequality also holds at the microscopic level, that is, for two different versions 
^l, S,i of the zero-range process, or even between Ut and fj. 

Our task will be to put this formal arguing into a rigorous proof. Before doing 
that, we need some definitions. We recall the formula for the norm in H-i(K) in 
terms of the Green function Q: for a function (or even a measure) u : K R, 

u{x)u{y)g{x, y)v{dx)v{dy). 

KxK 

More important for us will be the discrete version of this formula: for u : y„ ^ M 
such that u{ai) = 0, 



^ J2 u{x)u{y)g{x,y). 



M"'M-i,ri 22' 

x,y€Vn 

Now we define what we understand by "convergence in Let {j^"}n be a 

sequence of measures in ri„. Let u : K ^ [0, 00) be a given function. We say that 
I/" converges to u in the sense if 



lim / 

n—*oo J 



||e-u|P_i.„z^"(d?)-0. 
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A simple computation shows that the local equilibrium measures i^^(-) defined as 
the product measures in ri„ with marginals 

converge to u in the sense. 

For two given measures v, v' in , we define the relative entropy of v with 
respect to u' by 

I +CX), otherwise. 

We also say that v is stochastically dominated by v' if there is a measure A in 
il„ X ri„ such that 

i) A(e,a.) = i'(e) 

ii) \{n^,s) = v'{s,) 

iii) A{(C,r);e(a;) < for any a; G K?} - 1. 

Now we are ready to state the main result of this article. 

Theorem 4.1. Let he a sequence oj probability measures in , converging 

in the H^i sense to some bounded function uq : K [0, oo). Assume that the 
hydrodynamic equation has a unique solution. Assume also the technical conditions: 

i) Under (SG), there are positive constants k, p such that H{v^\i'p) < k3". 

ii) Under (C), there are two positive constants p < p' such that Vp is stochas- 
tically dominated by and i^" is stochastically dominated by Vpi for any 
n>0. 

Then, for any t > 0, the distributions {i^"(t)}n o-t time t > of the rescaled 
process ^" — ^s^t in f2„ starting from , converge in the sense to u{t,-), 

solution of the hydrodynamic equation 1^ . 

Under Condition ii), we have H{i'"\i'p) < H{vP \vp) and therefore there is a 
constant k such that H{v^\vp) < k3" for any n. 

4.3. The i7_i-norm method: martingale representation. In this section we 
will obtain a martingale representation for the _ff_i-norm of Notice that, due 
to the time-scaling, the generator of ^" is equal to 5"Lzr- For any function F : 
[0, T] X rifj*" M, differentiable in time and linearly growing in ^, Dynkin's formula 
states that 

Jo 

is a martingale. A long and tedious, but totally elementary computation shows 
that, in fact, 

where G{x) — Q(x,x) and 

H^,u)^{g{0~cly{u)){^-u)-g{0. 
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Observe that J !F{^{x),u)i'p{d(,) > for any p E [0, p*) and any m > 0. In 
particular, j T{£^{x),v)i''^^ -^{d£) > for any x G Vn- Therefore, 

is a martingale. Since Mq — 0, we have that E„[Af"] = for any t > 0. Here 
and below, P„ denotes the distribution of the process ^" starting from ly", and E„ 
denotes expectation with respect to P„. 

Taking the expectation with respect to P„ of the previous identity, we can obtain 
an expression for lE„||^t ~ """ll-i.n- 

E„||6~wril-i.„ = E«llCo-<|P-i,„-E„ f ^Yl ^iCi^), u^ix))ds 

/ i E 9{Ux))^nG{x)ds. 



Theorem |4.1| is an immediate consequence of the following two lemmas: 
Lemma 4.2. 



Lemma 4.3. 



lim En f ^{Cix),<\x))ds > 0. 

— Jo 

ft I 

lim E„ / ^ g{Ux))Ang{x)ds < 0. 



In fact, from these two lemmas, we conclude that 

lim sup E„ 1 1 - u J' 1 1 2_ ^ < lim sup E„ 1 1 Co - uo 1 1 ^ 1 , „ , 

n n — >oo 

and convergence in the sense follows at once. The argument behind the proof 
of these two lemmas is as follows. We will see that some sort of weak conser- 
vation of local equilibrium will allow us to replace in the previous expressions 
the functions g{^s{x)) by 0(5^ (a;)), where the symbol ^^(a:) denotes the average 
of S,siy) over a small triangle containing x, paying a price that vanishes when 
n ^ oo and then k oo. In the same way we can substitute !F{^s{x),Uf{x)) by 
(/)(Cs (x)) — (/)(m"(x))(Cs (x) — This last term is always positive, and we are 

in position to prove Lemma [4.2[ The proof of Lemma [4.3| is more subtle. Notice 
the huge factor 1/3" in front of the average in Lemma [4. 3 [ Since g{x,y) satisfies 
supj^ < -l-oo, we could guess that is in Hi{K). In that case, we 

should have A„5(x)/3" — > as n ^ oo in some convenient sense. It turns out that 
this is not the case. The function Q{x) is extremely irregular, and in fact it can 
be proved that for each x € V* fixed, A„5(a;)/3" 3/7 as n — > cx). Replacing 
g{^s{x)) by (/)(Cg (a;)) we will be able to perform an integration by parts in a small 
triangle of size k, therefore gaining a factor 3*^ that will save the day at the end. 

We will devote the following two sections to the proof of each one of these two 
lemmas. 
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5. The one-block estimate 

The replacement mentioned in the previous section is known in the hterature of 
interacting particle systems as the one-block estimate. Before stating the one-block 
estimate in a precise way, we need some definitions. 

Fix two integers n > ^ > 0. For x G Vn\Vn-i, we define Tj^{x) as the set of points 
in Vn contained in the triangle A(a;o, xi, 2:2), which contains x and has vertices in 
Vn-i- In an equivalent way, Tj^{x) = T OVn, where T is the unique triangle of the 
form (pi^_i o • • • o ipi^[K) containing x. For x S Vn-i, there are two possible choices 
for Tj^{x), given by two triangles intersecting exactly at x. Rotating the graph in 
such a way that both triangles lie on the upper semiplane, we choose T^{x) as the 
triangle at the right of x. The exact choice in this case is not important, the point 
is to choose each triangle the same number of times. 

Theorem 5.1 (Local one-block). Let us define the average number of particles 
isix) by 



where \Vk\ denotes the cardinality ofVk (and also ofTj^{x)). Then, 



lim limsup sup E„ 

^ n — >oo x£Vn 



ds 



= 0, 



(5.1) 



lim limsup sup 

k—^OD Yl^OO XGVn 



{U^)g{U^))-^{^^{x)){l+^':ix))}ds\=0. (5.2) 



The word "local" comes from the fact that in the usual version of the one- 
block estimate, the arguments in the integrals are averaging against a smooth test 
function. This local version of the one-block estimate was introduced in [5]. As 
stated in |5], this local one-block estimate is available only in dimension d < 2. 
The Sierpinski gasket K has Hausdorff dimension dn — log(3/2) < 2. But this 
is not really the point. The local one-block estimate holds each time the scaling 
of the process is faster than the scaling of the number of points. In our present 
situation, the process scales like 5" and the number of points scales like 3", so the 
local one-block estimate will hold. 



5.1. Proof of the one-block estimate. We will take the proof of Theorem |5.1| 
from [5j. In that paper the case on which condition (SG) is satisfied is treated in 
detail, so here we focus on condition (C). Our first step is to introduce a cut-off 
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function that prevents us to have too many particles at site x. Take a > p' . Then, 



)ds 



< 



{9iU^))-<f>{^':{x))}l{^^ix)>a)c 

< E„£ {giU^)) + cl,{e,{x))}l{es{x) > a)ds 

< E,, J* {g{U^)) + ^{C^{x))}l{^^{x) > a)ds 
<t f {g{i{x)) + <ly{e{x))}l{e{x)>a),y,,{dO 



< t J {g{^{x)r + 0(e'=(a:))'}v(rfOv(^'(^) > a)- 

The expectation in the last line is bounded in k. Moreover, by tha law of large 
numbers, £,''{x) converges to p' in probability as fc — > oo, and the probability in 
the last line goes to as /c — s- cx). Notice that this convergence is uniform in x. 
Therefore, we can introduce the indicatior function < a) in (5.1 1. 

By assumption, the entropy density _ff (j/"|i/p)/3" is uniformly bounded in n, 
by a constant k < +oo. A simple computation shows that the same is true for 
H^v^li^h), the entropy with respect to the invariant measure of the process. It is 
well known that the entropy of a jump process with respect to the invariant measure 
is decreasing in time. Fix some reference time T > 0. Denote by Pinv the law of 
the process up to time T, speeded up by 5" and starting from the invariant 
measure Vh- Then, there is another constant R depending only on k and T, such 
that H{Vn\Pinv) /i" < K for any n > 0. To simplify the notation, let us define 
Vk{^,x) by 

Vfe(e,x) = - ^{e{x))}lU^ix) < a). 

We have ommited in the notation the dependence of in n and a. By the 
entropy inequality, 

rT 



V{^s, x)ds 



< 



K 

< - 



7 78" 



logE, 



exp ■ 



{73"! f V,i^,,x)ds\j 



Using the elementary inequality e'^' < e 



we can get rid of the modulus 



in the previous expression. Therefore, the limit in (5.1 ) will be obtained if we prove 
that 



lim lim sup log E^. 

fe^oo „^oo 78" 



exp 



{±73"^ Vki^,,x)ds} 



0. 



By Feynman-Kac's formula, the logarithm of this expectation is bounded by the 
largest eigenvalue of the operator 5"'L ± 73"Vfe, where the term Vk is understood 
as a multiplication operator. For simplicity, we will consider just the "+" sign in 
Vk- By the variational formula for the largest eigenvalue of an operator in L'^{i'h), 
the previous expression is bounded by 



Tsup{(Vfc,/)-^(^)"(y7,-L,,y7)}, 



(5.3) 
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where the supremum is over all the densities / with respect to Vh, and the inner 
product is with respect to as well. Notice first that Vk depends on ^(y) only 
through its values for y e T^{x). Taking the conditional expectation of / with 
respect to !F{T^{x)), the a-algebra generated by {£,{y);y & T^{x)}, and due to 
the fact that i^h is a product measure, we can restrict the previous supremum 

to densities in the configuration space Nq" , which is homeomorphic to Ng*. By 
positivity of ( v7> ~LzrV7) i can also change L^r by the generator of a zero-range 
process restricted to the triangle T^{x). We will denote this generator simply by 
L, since no risk of confusion will appear by the fact that L depends on n, k and x. 
In this way, we have reduced the initial problem into a problem on a finite graph, 
which in our case is equal to Vk- Moreover, due to the presence of the indicator 
function l(^'^(a;) < a) in the definition of Vk, we can restrict ourselves to a finite 
state space, namely € Nq'';^'' < a}. Notice that at this point, ^'^ = ^^{x) does 
not depend on a: G 14. Since now the supremum is over a compact set, we can 
exchange the supremum and the limit as n — + cx3 to obtain 

limsupsupf(Vfc,/)--(^)"(v/7,-Lv/7>}= sup (V^,/). 

Now we have to identify the densities for which (v7j ~L\/f) = 0. A simple 
computation shows that, in this case, / is constant over the sets 

The restriction < a imposes I < a|Vfe|. Let us define the measures I'kj by 
taking Vk,i{') — Vp{'\^^ = ^l^fcD- Notice that these measures do not depend on the 
value of p, and they are also exchangeable. Then, the previous supremum is equal 
to 



sup / Vk{i,x)vk,i{dC)- 

l<a\Vk\ J 

For I < a|Vfc|, the indicator function < a) is identically equal to 1. There- 
fore, we are left with 

sup / {g{ax))-<P{l/\Vk\)}di^k,i. 

l<a\Vk\ J 

But this last quantity goes to as fc ^ cxd by the equivalence of ensembles, which 
states that the i^rancanonica/ measures i>k,i approach the canonica/ measures i^i/\Vk\j 
uniformly in compact sets of the real line. 

The case on which we take the "— " sign in front of Vk is totally analogous. In 



this way we have finished the proof of (5.1 1. The proof of this theorem when we take 



£.ix)g{^{x)) instead of 5(^(2;)) is entirely analogous, and we left it to the interested 
reader. 

Remark 5.2. The same estimate remains true if we consider functions of the form 
g(^"(x))F(t), where : [0,T] — > K is bounded. In that case the limit is uniform 
over sets of the form {supjgjQ |f (i)| < K}. It is enough to replace the variational 

formula in 



'5.S) by Xn{t)dt, where A„(i) is the largest eigenvalue of the operator 



5"L±73"Va: 
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5.2. Proofs of Lemma |4.2| and Lemma |4.3[ Now we are in position to prove 
Lemma |4.2[ By the one-block estimate, we can write 







plus a rest that vanishes as n ^ oo and then k ^ oo. But this last term is always 
positive, which proves Lemm a |4.2[ 

In order to prove Lemma 4.3 we start proving that A„t/(a;)/3" is uniformly 
bounded. Remember that Q(x,x) > Q{x,y) for any x,y Cz K. Since A„CJ(a;,a;) — 
—3", we conclude that Q{x^ x) < G{x, y) + 3" for x, y in Vn with y ~„ x. Therefore, 

Q{x,y) < G{x,x) < g{x,y), 
Q{x,y) < g{y,y) < G{x,y), 



where we have obtained the second line by interchanging the roles of x and y. We 
conclude that \G{x) — G{y)\ < 3" for x ~„ y, and |A„CJ(x)/3"| < 4. Now we are 
able to use the one-block estimate to rewrite the expectation in Lemma |4.3| as 

IE« E 9{U^))^ng{x)ds^E,, f 0fc'(a:))A„e(x)/3"ds 

■.'0 ^-1/ Jo 'J ^T/ 

plus a rest that vanishes as n ^ oo and then fc ^ oo. Notice that the function 
(j){^g{x)) is constant in T^{x). Therefore, we can integrate by parts (in this discrete 
context just a summation) the function A„5(x) in T^{x) to obtain that the previous 
expression is equal to 

En E '^fc'(^)) E dlkQ{^)ds, 

•^0 ' x6V„ 1=0,1,2 

where the symbols 9*^ kG{x) denote the outer normal derivative of Q{x) computed 
on the three vertices a"''^(a;) of the triangle T^^{x), defined by 

9;,,G(x) = (5/3)" i^^y) - 

y~n'\'''(.x) 

Notice that the ordering of the three vertices a^''^{x) is not relevant here. We 
have also gained a factor |Vfc|^^ in this integral. Now we just need to prove that 
the normal derivatives of Q {x) defined in this way are uniformly bounded. But the 
same arguments used to bound A„C7(a;) can be repeated here, to get a bound of 
the form kG{x)\ < 2 for any x £ Vn-k and any k < n. 

Appendix A. The behavior of the Green function at the diagonal 

In this Appendix we study the behavior of the function G{x). We want to 
compute A„t/(x)/3" for x £ Vn and we want to obtain its asymptotic behavior. 
The idea is to obtain an iterative formula for An^iQ(x) in terms of the values 
of G{x,y) for sites y in the neighborhood of x. In Fig. 1 we have taken a point 
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xo & Vn and we have named Xi, x'^, i = 1,2 each of the four neighbors of xq in 
Vn- Wc have also drawn points jji, y[, i = 0, 1,2, which correspond to the points 
in Vn+i belonging to the two triangles of side 2~" in Vn, meeting at xq- We claim 
that G{yi,Xj) can be computed in terms of the six numbers G^- = G{xi,Xj) (there 
are 9 combinations for /. but remeber that GJl- = G?,-)- In fact, let us construct 
G{y,yo) for y G Vn+1, for example. Defining 



Go{y,yo) 



3/10, y = yo 
1/10, y = yi,y2 
0, otherwise, 



we see that A„+i{3/5)"+^Go{y,yo) = -3"(5(t/, j/o), except for y = xi, i = 0,1,2. 
Therefore, the true Green function ^(y, j/o) is a linear combination between Go{y, yo) 
and G{y,Xi), i = 0, 1,2. In fact, A„+i(3/5)"+ieo(a;i, yo) = 3"+i/i°(a:i), where 



h°{x,) 



i = 
i = 1,2. 

Therefore, we have the following formula for G{y,yo)- 

.5. 



G{y,yo) = (ly^^Go{y,yo)+ J2 h\xi)G{y,Xi). 



i=0,l,2 



Notice, as well, that G{yi,Xj) can be computed from by harmonic continua- 
tion: 



Giy„x,) = l{2 J2 Gl,-G-}. 



fc=0,l,2 

Combining these two formulas, we can obtain the values of G{yi,yj) in terms of 
the numbers {G^} for any We give the formulas for G{yo,yo), G{yi,yo)\ the 
other formulas can be obtained by cyclic permutations of {0, 1, 2}: 

3 /3\"+i 1 

Giyo, yo) = (5) + ^{Goo + 4G1\ + 4G^2 + 4G^i + SG?^ + 4G^o}, 
G{yo,yi) = 1^ ( 5) +2^ {2G^o + 2Gri + 4G^2 + 5G^i + 6G?2 + 6G^o}- 
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Of course, similar formulas hold for the left-side points in Figure |A] As an appli- 
cation of these formulas, we can try to compute An^iQ(xQ) in terms of AnQ^xg): 

5"+^ {g{yi,yi) + Q{y2,V2) ~ G[?„) = 3"+^ • \+ 

+ 5"-i{5(Gl\ + G^i^ - 2GJf„) + 12(G;,\ + GJJ^ - G«o) + 8(G?2 - Goo)}. (A.l) 

Adding the symmetric term coming from the left-hand side of Figure |X] we see 
that 

An+i^K) = 3"+i • I + A„^(xo) + ^A„g(xo,xo) 

+ 5" • \{Q{x^.X2) + Q{x\,x'^) ~ 2g{xQ,Xo)]. 
In this expression, there is a new term appearing. Let us define 
r„(a;o) = 5"{c;(a;i,a;2) + Q{x\,x'2) - 2Q{xq„Xq)). 
Since A„5(a;o,xo) = —3", we have the formula 

An+iQ{xa) = 3"+i • ^ + A„^(a;o) + ^r„(xo). 

Now let us compute r„+i(a;o): 

r„+i(a;o) = -3"+i • i + ^A„g(a;o) +r„(a;o). 

This establish a linear recurrence formula for the pair {A„t/(a;o), r„(a;o)}- Due 
to the factor 3"^^, it is natural to define a„ = A„ (a;o)/3". For 

a„, bn, the recurrence formula reads 

a„+i = 2/5 -I- a„/3 + 8&„/15 
bn+i = -1/5 + 2a„/15 + 6„/3. 



This recursion formula can be written in vectorial terms as a„+i = w + Ma„, 
with 

M-( ^/^^ 
" 2/15 1/3 

The eigenvalues of this matrix are Ai — 3/5, A2 — 1/15. In particular, for any 
initial value of a„, &„ (remember that a;o G Vn, so the sequence does not start at 
n = 1) we have convergence to a unique fixed point, given by a = (/ — il/)^^w. In 
our case, w= (2/5,-1/5), and a = (3/7,-3/14). In particular, 

lim A„5(2;o)/3" = 3/7. 

n — *oo 

Notice that the limit is positive and does not depend on xq. Remember that 
A„t/(a;o,a;o)/3" = —1. This remarkable fact shows the high irregularity of the 
function G{x). In contrast to it, for the unit interval [0, 1], G{x) = x{l — x), so G{x) 
is smooth and concave. 
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In order to obtain a recursive formula for the partial derivatives of Q{x), we focus 
our attention on the right side of Figure |A] Let us take a look at equation ( A.l I. 
We see that most of the work has already been done. In fact, defining 

a„ = (5/3r(G?i + G^2-2G^y 

7„ = (5/3r(GSi + G^o-2G^o), 
we obtain the following recursion formulas: 

a„+i = 3/5 + l/3a„ + 8/15/3„ + 4/57„ 
(3n+i - 1/10 + 2/15a„ + l/3/3„ + 2/57„ 
7„+i = 7„. 

In particular, we obtain the same linear recursion as before, but now with a 
difi^erent vector w = (3/5 + 4/57, 1/10 + 2/57). Again, we have convergence of 
(a„,/?„) to the solution of a = w + Ma. In our case, lim„ a„ = a, with a — 
17/14 + 27, where 7 = 7„ for some suitable n. Notice that we recover our previous 
computation A„t/(a;)/3" — > 3/7 by noticing that the Laplacian at x is the sum of 
the two partial derivatives at x, and that the corresponding 7's add up to —1 in 
that case. 

We see that the behavior of AnG{x) is the worst possible, in the sense that in 
one hand we have the trivial bound A„CJ(a;)/3" < Ci for any n and Ci = 4, and in 
the other hand we have the lower bound A„tJ(x)/3" > C2 for any C2 < 3/7 and any 
n large enough. 
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